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ABSTRACT 

Damped transverse oscillations of magnetic loops are routinely observed in the solar corona. This phenomenon is interpreted as 
standing kink magnetohydrodynamic waves, which are damped by resonant absorption owing to plasma inhomogeneity across the 
magnetic field. The periods and damping times of these oscillations can be used to probe the physical conditions of the coronal 
medium. Some observations suggest that interaction between neighboring oscillating loops in an active region may be important 
and can modify the properties of the oscillations compared to those of an isolated loop. Here we theoretically investigate resonantly 
damped transverse oscillations of interacting non-uniform coronal loops. We provide a semi-analytic method, based on the T-matrix 
theory of scattering, to compute the frequencies and damping rates of collective oscillations of an arbitrary configuration of parallel 
cylindrical loops. The effect of resonant damping is included in the T-matrix scheme in the thin boundary approximation. Analytic 
and numerical results in the specific case of two interacting loops are given as an application. 
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1. INTRODUCTION 

Transverse oscillations of magnetic loops in the solar corona 
are under intense research since the first observational reports 
by th e Transition Region And Coronal Explorer ( TRAC E) (see, 
e.g., iNakariakov et al.1 Il999t lAschwanden et all Il999h . Large- 
amplitude coronal loop oscillations are usually excited af¬ 
ter energetic events as, e.g., solar flares, coronal mass ejec¬ 
tions, or low coronal eruptions (see IZimovets & Nakariakovl 
20151) . Based on magnetohydrodyna mic (MHD) wave theory 
(e.g., INakariakov & Verwichtel l2005h . transverse loop oscilla¬ 
tions have been interpreted as standing kink MHD waves. 
Kink MHD modes are nearly incompressible waves, mainly 
driven by magnetic tension, and responsible for global trans¬ 
verse motions of the flux tube (see, e.g., lEdwin & Roberts 
1 9831 iGoossens et al.ll2009[ 120121) . A relevant feature of large- 
amplitude loop oscillations is that they are strongly damped. It 
has been shown that resonant absorption, caused by plasma in¬ 
homogeneity in the direction perpendicular to the magnetic field, 
is an efficient damp ing mechanism of kink MHD waves in coro¬ 
nal lo ops (see, e.g.. lRuderman & R oberts 2 0021: IGoossens et~aU 
120021) . Due to resonant absorption, the energy from the global 
kink motion of the flux tube is transferred to small-scale, unre¬ 
solved rotational m otions around the nonuniform boundary of 
the tube (see, e.g., Teiradas et al.1 l2006t IGoossens et al.1l2014t 
ISoler & Te rradas 2015). As a result of this process, the global 
kink oscillation of the coronal loop is quic kly da mped in time. 
The interested reader is referred to lGoossens et al . (1201 ll) . where 
the theory and applications of resonant waves in the solar atmo¬ 
sphere are reviewed. 

Observations often show that neighboring oscillating loops 
in an active region interact with each other and exhibit collective 


behaviour (e.g., Schriiver & Brown 

2000 

; [Verwichte et alJ2004; 

Schriiver et al. 2002j White et al. 

2013 

). Interaction between 


loops can modify the properties of their transverse oscillations 
compared to those of the classic kink mode of an isolated loop. 
Therefore, advanced models describing coronal loop oscillations 
should take into account interactions within loop systems. A 
number of works have studie d collective transverse oscillations 
in Cartesian geometry (e.g.. IDfa z et al.1 20051: iDfaz & Rohertsl 


120061: iLuna etal. 2006|: Arregui et akf 20071 l2008h . In cvlindri 


cal geometry, Luna et al.l (120081) numerically inv estigated trans 


verse oscillations of two cylindrical loops, and lOfmanl {2009) 
performed numerical simulations in the case of four interact¬ 
ing loops^ Concerning analytical works in cylindrical geome¬ 
try, iLuna et aU {2009L 2010|) us ed the T-matrix theory of scat¬ 
tering (see. e.g.. lTwerskvl 19521: IWatermanll 19691: lBogdan|[l987l: 


Ken nens et al.lll994l) to investigat e transverse oscillations of two 

and three loons (ILuna et al .1 12009 ) and of bundles of many loops 
(ILuna et al.ll2010l) in the B = 0 approximation, where /? refers to 
the rat io of the gas pressure to the mag netic pressu re. ISoJereLalJ 
(120091) later extended the method of iLunai et al.l (120091 2010) 
by incorporating gas pressure and longitudinal flows and stud¬ 
ied collective oscillations of flowing prominence threads. These 
works showed that loop interaction affects the properties of their 
oscillations. ILuna et all ( 2008 , 2009) obtained that a system of 
two loops of arbitrary radii supports four kink-like collective 
modes. The shift of the collective mode frequencies with respect 
to the individual kink frequencies of the loops is significant when 
the distance between loops is small (of the order of the loop ra¬ 
dius) and when loops have similar densities. Conversely, the os¬ 
cillating loops show little interaction when they are far from each 
other and when their densities are substantially different. On the 
other hand, IVan Doorsselaere et al.l ( 2008 ) and iRobertson et al.l 
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( 20id ) used a different method based on bycilindrical coordi¬ 
nates to study transverse oscillations of two pressure-less loops 
in the thin tube (TT) approx imat ion. Of the four kink-like collec¬ 
tive modes obtained bv ILuna et al. ( 2009 ) in the T-matrix theory, 
only two different modes remain in the TT approximation con¬ 
sidered by IVan Doorsselaere et alJ (i2008i) and iRobertson et al.l 
(1201 Oh . 

Concern ing the damping of the oscillations, lArreeui et alJ 
(2007, 2008) in vestigated resonan tly damped oscillations of two 
slabs, while iTbrradas et akl (|2008) numerically studied the res¬ 


onant damping of transverse oscillations of a multi-stranded 
loop. Those works showed that the process of resonant damp¬ 
ing is not compromised by the irregular geometry of a real¬ 
istic loop model and still produces the efficient attenuation of 
global transverse oscillations. The damping of transverse oscil¬ 
lations of two cylindrica l loops was analytically investigated by 
IRobertson & Rudermanl d201 If) and iGiisen & Van Doorsselaere! 
(120141) . who considered the TT approximation and used bycilin¬ 
drical coordinates. Results obtained with bycilindrical coordi¬ 
nates should be treated with caution when the distance between 
loops is small. Geometrical effects intrinsically associated to 
the bycilindrical coordinates may produce unphysical results. 
Our purpose is to use the T-matrix method of iLuna et alJ (120091 
1 2010 ) to investigate resonantly damped oscillations of bundles 
of loops. The present paper impartially based on unpublished 
results included in Solen (1201 OtFl The effect of resonant ab- 
sortion in the Alfven continuum is incorporated to the T-matrix 
scheme by using the method that combines the jump conditions 
of the perturbations at the resonance position with the so-called 
thin boundary (TB) approximation (see, e.g. JSakurai et alll99lb 


iGoossens et al. 


_app 

119921). A similar method has previously been 
used by iKeppensl ( i 1995b to investigate absorption of acoustic 
waves. We provide a general analytic theory, which is valid for 
bundles of many transversely nonuniform parallel loops of ar¬ 
bitrary radii. Specific results in the case of two loops are ob¬ 
tained and compared to those given in [Robertson & Rudermanl 
(1201 ll) and lGiisen & Van Doo rsselaere] (120141). Our r esults are 
also compared to those of lArregui et al.l (120071120081) obtained 
in Cartesian geometry. 

This paper is organized as follows. Section[2]contains the de¬ 
scription of the equilibrium configuration and the basic govern¬ 
ing equations. The general analytic T-matrix theory of scatter¬ 
ing to compute the frequencies and damping rates of collective 
loop oscillations is given in Section[3] Later, the specific case of 
damped oscillations of two loops is discussed both analytically 
and numerically in Section [4] Finally, some concluding remarks 
are given in Section 0 


2. MODEL AND GOVERNING EQUATIONS 

Our equilibrium configuration is composed of N straight and par¬ 
allel magnetic cylinders of length L embedded in a uniform coro¬ 
nal plasma. The ends of the magnetic tubes are fixed at two rigid 
walls representing the solar photosphere. We set the ^-direction 
to be along the axes of the tubes. The magnetic field is straight 
along the z-direction, namely B = Be z , where B is a constant ev¬ 
erywhere. We use subscripts ‘i’ and ‘e’ to refer to, in general, the 
internal region of the tubes and the external plasma, respectively. 
The subscript or superscript ‘j’ is used to refer to a particular 
loop. We denote by R i the radius of the jth tube. The distance be¬ 
tween the centers of the jth and j’th loops is d^. We denote by pj 



Fig. 1. Sketch of the equilibrium configuration in the specific case of 
two transversely non-uniform coronal loops. 


the internal density of the jth tube, while p e denotes the external 
density, i.e., the density of the coronal environment. In our model 
pj and p e are constants. There is a transversely nonuniform tran¬ 
sitional layer surrounding each magnetic tube in which the den¬ 
sity continuously varies from the internal density, pj, to the ex¬ 
ternal density, p e . The thicknesses of the non-uniform boundary 
layer of the jth cylinter is /j. A sketch of the equilibrium config¬ 
uration in the case of two magnetic tubes (N = 2) is given in 
Figure Q] 

We adopt the (5 — 0 approximation, where /) refers to the ra¬ 
tio of the thermal pressure to the magnetic pressure. This is an 
appropriate approximation to investigate transverse waves in the 
solar corona. In the (3 = 0 approximation, the ideal MHD equa¬ 
tions governing linear perturbations superimposed on the static 
equilibrium state are 



- (V x b) x B, 
B 

V x (£ x B), 


(1) 

( 2 ) 


where £ is the plasma Lagrangian displacement, b is the mag¬ 
netic field Eulerian perturbation, p is the density, and p is the 
magnetic permittivity. 

We assume the temporal dependence of perturbations as 
exp {—iiot), where to is the oscillation frequency. In the case of 
transversely nonuniform tubes, the global transverse oscillations 
are quasi-modes whose frequency is complex owing to damping 
by resonant absorption, i.e., to = wr + ion, where wr and to\ 
are the real and imaginary parts of the frequency, respectively. 
The real part of to is related to the period and the imaginary part 
corresponds to the damping rate of the oscillations. We consider 
that the oscillating flux tubes are line-tied at the photosphere, 
which acts as a perfectly reflecting wall in this model owing to 
its large density compared to the coronal density. Hence, we as¬ 
sume perturbations to be proportional to exp (ik z z), with k the 
longitudinal wavenumber. For standing oscillations, k z given by 

1171 

k z =—, with n=l,2,... (3) 

L 


1 The full text of ISoled d2010h is available at We shall restict ourselves to the fundamental mode of oscillation, 
http://www.uib.es/depart/dfs/Solar/thesis_robertosoler.p<Mj we take n = 1. 
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In the regions with constant density. Equations (T} and ([2} 
can be reduced to following equation, 

vi p' + kip' = o, ( 4 ) 


where P' = B • b//r is the total pressure Eulerian pertubation 
and the subscript ± refers to the direction perpendicular to the 
magnetic field. Thus, Vy denotes the perpendicular part of the 
V 2 operator. In turn, the quantity k± plays the role of the perpen¬ 
dicular wavenumber and is defined as 


k 


2 

_L 



(5) 


where co 2 A = k 2 v 2 is the square of the Alfven frequency and 
= B 1 1 up is the square of the Alfven velocity. We stress that 
Equation © is only valid in the regions with constant density, so 
it does not apply within the nonuniform boundaries of the loops. 


3. T-MATRIX THEORY OF SCATTERING 


Equation Q is the two-dimensional Helmholtz Equation. To 
solve Equation (4]) we use the scattering theory in its T-matrix 
formalism (see, e. g.. 1 Watermanl [19691) . In the solar context, the 
T-matrix theory has previously been used to investigate the scat¬ 
tering and absorption properties of bundles of magnetic flux 
tubes (e.g. Bogdan & Zweibelll 19851: iBogdan & Cattaneol 19891: 
iKeppemetal. 11994 : Kepnenslll995 , among others PLuna et al l 
( 20091 120 1 0t> used of this technique to compute the eigenmodes 
of systems of magnetic tubes. Because of the inhomogeneity of 
the tubes in the transverse direction, the modes with frequen¬ 
cies between the internal Alfven frequencies of the loops and 
the external Alfven frequency are resonant in the Alfven con¬ 
tinuum. As a result, the oscillations are damped by resonant ab¬ 
sorption. The effect of resonant absorption was not considered 
bv lLuna et alJ ( 2009 , 2010 !). Here we extend their theory to con¬ 
sider resonant damping. 


3. 1. Solutions in the Internal and External Plasmas 

We use local polar coordinates associated to the jth loop. We 
denote by r t and cp t the radial and azimuthal coordinates, respec¬ 
tively, of the coordinate system whose origin is located at the 
center of the jth tube. We can define an equivalent coordinate 
system in each tube. In this local coordinate system, the solu¬ 
tion to Equation (J4]) in the internal region of the jth tube can be 
expressed as 

oo 

P'i = 2 A\„j m (^xpi) exp (inup^j , (6) 

m——oo 


waves. The total external solution is computed by adding the net 
contributions of all flux tubes, namely 

^ = 2>e j , (7) 

j 

where P'l is the net contribution of the jth tube to the external 
solution. The key idea behind the scattering theory is that the so¬ 
lution of Equation (|4ji in the external plasma can be decomposed 
into several fields with different physical meanings, namely the 
total , exciting , and scattered fields. Here we give an overview of 
the method. Interested readers are referred to [Luna et alJ (20091. 
12010 ) for extensive explanations. 

In the external plasma, the total field associated to the jth 
cylinder can be expressed as 

OO 

P total= 2 + °^2,n H m { k ^ r t)] eX P (//»^j) , (8) 

m——oo 

where //)„' 1 and H i2> are the usual Hankel functions of the first 
and second kind, respectively, and a\ m and tr\ ; constants. The 
first term on the right-hand side of Equation © represents out¬ 
going waves from the jth tube, and the second term represents 
incoming waves toward the jth tube. Importantly, we note that 

K*Y, P 'Lr (9) 

i 

The reason for this inequality is that the outgoing wave asso¬ 
ciated to a particular tube contributes as an incoming wave for 
all the other tubes. In other words, £j P[\ )lA is not the net exter¬ 
nal solution. To overcome this problem, the total field associated 
to the jth cylinder, ,, is decomposed into a scattered field , 

/’scat’ and a exciting field. P'* xc ; . Conceptually, the scattered field 
represents the actual contribution of the various tubes to the net 
external solution, whereas the exciting field can be understood 
as the cross-talk mechanism r esponsible for interaction b etween 
flux tubes (see details in, e.g jBogdan & Cattanedfl989l) . 

The full net external solution, P' e , is defined so that it corre¬ 
sponds to the sum of the scattered fields associated to all tubes, 
namely 

^ = 2>slar (10) 

j 

Conversely, the exciting field associated to the jth tube is defined 
as the difference between the full net contribution and the scat¬ 
tered field of the jth tube, namely 

p'l _ p' _p'J -V pm (-JJ-v 

excit r e ’scat ^/j scar 
n#j 


where m is the azimuthal wavenumber, J m is the usual Bessel 
function of the first kind of order m, and A],, are constants. Unlike 
the case of isolated tubes (see, e.g jEdwin & Roberts! 19831) , the 
solution is not entirely described by a single value of in. Because 
of interaction between tubes, the values of m are coupled. For 
transverse, kink-like oscillations the dominant terms in the ex¬ 
pansion are those with m = +1, but the contribution from other 
nisjs not negligible unless the tubes are far from each other 
(Luna et al.ll2009l) . 

The solution to Equation (4]) in the external region is written 
using the principle of superposition, which is applicable to linear 


iWaterman! ( 1969 ) introduced the T-matrix operator, T J , which 
linearly relates the scattered and exciting fields as 


p'j — r r) p'j 

scat excit" 


( 12 ) 


Bogdan ( 19871) show ed that for cylindrica l scatterers the T- 
matrix is diagonal, and Keppens et al. (1994) gave an expression 
of its elements, namely 


T J — 

' mm 


“Li > 


(13) 
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where a\ m and a\ m are the same constants the appear in Equa¬ 
tion ©. We can use Equation (fl3l > to eliminate ar’j and write all 

the expressions in terms of tri, m alone. With the help of these last 
formulae, and after some algebraic manipulations using well- 
known properties of the Bessel functions, we can rewrite Equa¬ 
tion © as 

oo 

Cal = Z 2(r L\ J « (^c r i) - CC (^etj)l exp (imipj ), 

m=-o o 

(14) 

from where it is straightforward to identify both exciting and 
scattered fields, namely 

OO 

Cat = Z 2ct , 2m J m (k ±c r i )exp(im<Pi), (15) 

m=— oo 
oo 

Cat = - Z 2( C 7 ’CC ) (^xcfj) exp . (16) 

m——oo 

Finally, we use the expression of P' 3 c . dl into Equation ( flOb to arrive 
at the total net solution in the external plasma, namely 

OO 

= - Z Z 2< 4 m T «» nH m ) { k ^ r i) eX P • (17) 

j m——oo 

Equations © and ( fTTb formally describe the total pressure 
perturbation in the interior and in the exterior of the tubes, re¬ 
spectively. However, we recall that these expressions do not ap¬ 
ply in the nonuniform boundary layers. The T-matrix elements, 
T 3 mm , contain the information about how the solutions are con¬ 
nected across the non-uniform boundaries of the tubes. 


3.2. T-matrix Elements in the Thin Boundary Approximation 


At this stage we incorporate the effect of the nonuniform bound¬ 
ary layers. In the nonuniform layer of thejth tube the global wave 
modes are resonant in the Alfven continuum at the resonant po¬ 
sition, rj = r Aj , where the global oscillation frequency matches 
the local Alfven frequency. The resonant position, r A| j, is de¬ 
fined through the resonant condition or = &?v^(r A j). We use 
the TB approximation and restrict ourselves to lj/Rj « 1. The 
TB approximation assumes that the jump of the perturbations 
across the resonant layer is the same as their jump across the 
whole nonuniform layer. Thus, the connection formulae of the 
wave perturbations across the resonance are used as jump con¬ 
ditions for the total pressure and the Lagrangian displacement 
at the boundaries of the tubes. This method and its applications 
have been reviewed by iGoossens et ali ( 2011) . The TB approx¬ 
imation within the formalism of the T-m atri x theory has prev i- 
ously been used bv lKeppens et al.l (119941) an d lKeppensl (! 19951) . 

General expressions of the connection formulae for the 
perturbations across the resonant layer can be found in, e.g., 
ISakurai et alj (1199 ll) . In local coordinates, the connection formu¬ 
lae for the total pressure Eulerian pertubation, P', and the radial 
component of the Lagrangian displacement, , at r, = /- A | are 


m = o. 


&] = 


. m2/r li p , 

m - P 

lpA A |j 


( 18 ) 


where [A] - X e - Xj denotes the jump of the quantity X across 
the resonant layer and |pA a Ij is defined as 



where we have used the resonant condition to 2 — k:v\(r A] ). Af¬ 
ter imposing the jump conditions given in Equation (fl8l) . we ob¬ 
tain both the T-matrix elements and the dispersion relation of the 
collective modes. 

On the one hand, the T-matrix elements are 

kir JUX'R,) _ k ±l j;„(k ±i R,) mT rj, 

^ _ Pe(^-A-M.) J,„(k ±c R j) J m (k ±i R,) IpAaI 

■C(ftxAi) , ■ 

J '"( k ±i R i) ^ AaI 

( 20 ) 

where the prime ' denotes the derivative of the Bessel or Hankel 
function with respect to its argument. In the absence of resonant 
damping. Equation (l20l > consistently reverts to Equation (17) of 
iLuna et al . (2009). 

On the other hand, the constants a J 2m satisfy an algebraic sys¬ 
tem of equations, namely 

OO 

4m + Z Z (*re</jj') ex P (i(nT - m)ipff) = 0, 

j'^j m'=—oo 

( 21 ) 

for -oo < m < oo, where tpjy is the angle formed by the vec¬ 
tor positions of the centers of the two tubes with respect to the 
global reference frame. Equation (l2TT > is system with an infinite 
number of algebraic equations. The condition that there is a non¬ 
trivial solution, i.e., the determinant formed by the coefficients 
set equal to zero, provides us with the dispersion relation (see 
details in iLuna et al.l 120091. l 2010h . The solution of the disper¬ 
sion relation is the complex frequency of the damped collective 
quasi-mode. The imaginary part of the frequency is the reso¬ 
nant damping rate. In the case of ILuna et all ( 2009 . [2010 }). the 
solution frequency was real because of the absence of resonant 
damping. 

The dispersion relation is a complicated expression and has 
to be solved by numerical methods. For practical computational 
purposes, in Equation (l2TT > the indices m and in' must be trun¬ 
cated into a finite number. The truncation term must be large 
enough to avoid inaccuracy of the solutions. Then, the disper¬ 
sion relation can be solved using standard numerical routines to 
find the roots of transcendental equations. 


*„■ ^'"'(AcRi) 

Pei^-kh’U H^(k ±e Rj) 

Jm (krcRi) 
(*Le*j)’ 


4. APPLICATION TO DAMPED OSCILLATIONS OF 
TWO LOOPS 

We recall that the T-matrix theory described in the previous Sec¬ 
tion is valid for an arbitrary number of loops. Here we apply that 
method to investigate damped transverse oscillations in a loop 
system composed of two magnetic tubes alone. First we derive 
approximate expressions of the period and damping time in the 
case of two identical thin tubes. Later, we consider two tubes 
with different properties and arbitrary radii and perform a nu¬ 
merical study. 

4.1. Approximate solutions for two identical thin tubes 

We consider the paradigmatic case of two tubes with identical 
properties. We denote by p x , R, and / the internal density, the 
radius, and the thickness of the boundary layer of both tubes. 
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respectively. The external density remains p e . The distance be¬ 
tween the centers of the two tubes is d. 

To obtain an approximate dispersion relation for kink-like 
modes we assume that the contribution from the azimuthal 
wavenumbers m — +1 to the collective oscillation is much more 
important than the contributions from other values of m. Thus, 
we take m, ml = ± 1 in Equation d2TT> and roughly neglect the 
other terms. As pointed out bv lLuna et al.l d2009l) . the coupling 
between the different values of m gets stronger as the separation 
between the cylinders decreases. In terms of the parameters of 
the model, this means that our approximation applies to the case 
of large separations, i.e., R/d <sc 1. 

The system defined in Equation (l2H now becomes an al¬ 
gebraic system of four equations for the unknowns cr\ t , a\ p 
a\_ v and a 2 r We obtain the dispersion relation from the con¬ 
dition that there is a non-trivial solution of the system, namely 


4.1.1. Frequency of the Oscillations 

First we neglect the presence of the nonuniform boundary lay¬ 
ers and study undamped oscillations. We drop the last term on 
the left-hand side of Equation (l25l > and the dispersion relation 
simplifies to 


Pi [co 2 - k 2 z v\) 6+p e [or - k 2 v 2 Ac ) = 0. 


(27) 


The exact solution is 


rut 


1 *£(§)“ 


(28) 


where £ = pi/p e is the density contrast and coi js the frequency 
of the kink mode in an individual thin tube (lEdwin & Roberts! 
1983 ). namely 


(■ K l) (k ±e d ) 2 - H^(k ±e d) 2 f T 4 u 

- 2 ( H ( 0 l) (k ±e d ) 2 + H^\k ±e d) 2 ) T 2 n + 1 = 0, (22) 

where we used the property that T\\ = T_|_| according to 
the symmetry relations of the Bessel and Hankel functions (see 
lAbramowitz & Stegunlll972l) . Equation (l22l > is an equation for 
co, which is enclosed in the definitions of T\\ and k ±e . To go 
further analytically, we perform a series expansion of the Han¬ 
kel functions for small arguments and retain the first term alone. 
This would be approximately valid for thin tubes. After some 
algebraic manipulations. Equation (l22l > can be recast as 

T u ± i^(k ±e d) 2 * 0. (23) 


Now, we look for a simplified expression of T\\. We use the 
TT approximation, i.e., R/L « 1. In Equation < 12()b we perform 
an asymptotic expansion of the Bessel and Hankel functions for 
small arguments. Equation (l20l i becomes 


_1_!__ + _ nlE _ 

p e (or-kh’ 2 Ac ) w 2 |dp/dr| fi 

__!_ + _!_ ijr _ HLH _ 

Pe(w 2 -A-Ji^ e ) Pi^-kh’ll) a>l\Apldr\ R 

in (k ±e R) 2m 
2 2m m\(m - 1)! ’ 


(24) 


where we also assumed ~ R for simplicity. We take m = 
1, substitute Equation (l24l > into Equation < 1 27b . and arrive at the 
following expression. 


Pi'ii+P ' V ie ,2 

Pi + Pe Z 


(29) 


Equ ation (l28l > can be 


IVan Doorsselaere et alJ (2' 


compared to Equation (51) of 
. Both expressions agree if the pa¬ 


rameter E of IVarTDoorsselaere et al. (20081) is approximated as 


E = exp[-2arcco sh(^)] ~ , which i s valid for R/d « 1. 

In agreement with I Van Doorsselaere et al. (2008). in the TT ap¬ 
proximation we obtain two kink-like solutions corresponding to 
the + and - signs in Equation (l28l >. As explained bv lLuna et al.l 
( 2009 ). there are actually four kink-like modes beyond the TT 
approximation. For arbitrary radii, the low-frequency solution 
(that with the + sign in Eq uation (l28l >) splits in the S x and A y 
modes of iLuna et al.l (20091), while the high-frequency solution 
(that with the - sign in Equation (l28l >) becomes their S y and A x 
solutions. 

From Equation (l28l > we compute the period of the oscilla¬ 
tions, P = In/u), as 



(30) 


with Pk = 2n/oJt the period of the kink mode of an isolated 
tube. For large separations the tubes feel little interaction. Then, 
the periods of the two solutions consistently tend to that of the 
kink mode of an isolated tube. 


4.1.2. Damping Rate 


p; (w 2 - k 2 z v 2 M )6 + p e (w 2 - k 2 v i e ) 

. n PiPe (^ 2 ~ k z v li) {u 2 ~ frMe) 

l R\dp/dr\ R co 1 

where the parameter <5 is defined as 


(25) 



Equation (l25l > is the dispersion relation for resonantly 
damped kink modes of two identical tubes in the TT and TB ap¬ 
proximations and for large separations. In the limit that the tubes 
are far from each other, R/d—) 0 and so 6 —> 1. Then, Equa¬ 
tion (l25l > consistently reverts to the dispersion relation of reso¬ 
nant kink modes in an isolated thin tube (e.g., iGoossens et al.l 
1992). 


Here we incorporate the effect of resonant damping and consider 
the full expression of the dispersion relation (Equation (l25l>). To 
obtain an approximate expression of the damping rate, we sub¬ 
stitute co = + ioj\ in Equation (l25l) and assume weak damping, 

i.e., |cui| <sc |c<jr|. Then, we are allowed to neglect terms of order 
coj and higher orders. After some algebraic manipulations (not 
given here for the sake of simplicity) we arrive at an expression 
for the ratio co\/co^, namely 


Oh n 1 P,p e 1 

WR 2Rpi6 + p e \dp/dr\ R co\ 

To simplify Equation (PTT) we consider a smooth monotonic pro¬ 
file for the density in the nonuniform boundaries so that we can 
express the derivative of the density profile as 


= kk 2 pi-p z 

dr R 4 l 


(32) 
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with T a factor that depends on the form of the profile. For in¬ 
stance, F = A/n 2 for a linear profile and T = 2fn for a sinu¬ 
soidal profile. In addition, we approximate wr by the value of 
the frequency in the undamped case (Equation (l28l >) and substi¬ 
tute the expression of 6 (Equation (26l >) to arrive at 


wr 


1 1 l £-1 
~2jrT~Rt/ + 1 





(33) 


For a linear density profile. Equation (l33l i becomes Equa¬ 
tion (86) of lRobertson & Rudermanl (1201 ll) if the approximation 

exp[-2arccosh(^)] ~ valid for R/d <sc 1 is again per¬ 
formed in their expression. 

Now, we define the damping time as td = 1 /\oj\\ and use 
Equation (l33l > to obtain the expression for the ratio r n /P. By 
keeping terms up to (R/d) 2 , the expression for the damping ratio 
is 


IE, (IE) [j T JL(*) 2 

P \P/k[ £+1 \d) 


(34) 


where 


r R {+1 
\pk i { - r 


(35) 


is the damping ratio of the individual kink mo de (see, e.g., 
iRuderman & Robertsll2002t iGoossens et al.ll2002l) . As for indi¬ 
vidual kink modes, the damping ratio of collective modes is in¬ 
versely proportional to the thickness of the nonuniform layer. 
Equation (34l > shows that the solution with the + sign is more 
efficiently damped by resonant absorption than the solution with 
the - sign. This difference in the damping rates of the two rnodes 
qualitatively agrees with the results of lRobertson & Rudermanl 
( 2011 ). 


4.2. Numerical Study 

Here we perform a numerical study beyond the approximate lim¬ 
its studied before. To do so, we consider the full dispersion rela¬ 
tion obtained from Equation (fTit by using the general expression 
of the T-matrix elements (Equation (20l) ). The truncated disper¬ 
sion relation is then numerically solved. Convergence tests of 
the results have been performed to make sure that the truncation 
value of the azimuthal series is large enough for the error in the 
solutions to be negligible. In short, we found that for kink-like 
modes the truncation value, namely m t , has no important effect 
unless m, < 5 and the two loops are next to each other (d/R ~ 2). 
When in, > 5 and/or the loops are far from each other, we essen¬ 
tially obtain that the results are independent of m t . We have used 
m, = 30 in all computations given here, which reduces the error 
and assures the excellent converge of the solutions even when 
d/R ~ 2. 

The numerical method follows a two-step procedure. First, 
we solve the dispersion relation in the absence of nonuniform 
boundary layers. In that case, the solution is a real frequency 
that corresponds to an approximation to wr. This approximate 
value of wr is used to compute the resonance positions and the 
derivative of the density profile at the resonances. We assume 
sinusoidal density profiles in the nonuniform boundary layers. 
Then, we use these parameters to solve the complete dispersion 
relation, which now includes the effect of resonant damping. The 
frequency obtained from the second run is complex, so that it 
provides us with a more accurate value of wr and also gives us 
the value of ( 0 \. 


4.2.1. Identical loops 


We initially study the case of two identical tubes. We set the 
Cartesian coordinates system so that the xy-plane is perpendicu¬ 
lar to the axes of the loops. In that plane, the centers of the two 
loops are located on the x-axis. We use the notation introduced 
bv lLuna et alJ d2008l) to denote the four kink-like modes present 
in a two-loop configuration. The modes are labeled as S x , A x , 
S y , and A y , where S and A denote symmetric or anti-symmetric 
motions of the two loops, respectively, and the subscripts x and 
y indicate the main direction of polarization of the oscillations 
in the coordinates system defined above. The eigenfunctions of 
these four modes in the case of loops without nonuniform bound¬ 
ary layers can be found in Figure 2 of lLuna et alJ ( 2008 ). 

Figure [2] shows the dependence of wr and the ratio |wi| / 
on the separation between loops, d/R, for a particular set of 
parameters given in the caption of the figure. We note that the 
curves corresponding to the S x and A y modes, and those of the 
A x and S y modes, are almost superimposed because we are in the 
TT regime (we used L/R = 100). Concerning the behavior of wr, 
the frequencies of the four solutions tend to the kink frequency 
of an isolated loop in the limit d/R » 1. Conversely, the smaller 
the separation between loops, the more important the splitting 
of the collective frequencies with respect to the kink frequency 
of an isolated loop. Figure [2 a) can be compared to Figure 3 of 
iLuna et all (20081) and with Figure 4 of [Van Doorsselaere et al.l 
( 20081) . However, unlike in those previous works we note that in 
our case the frequencies of the high-frequency solutions (A x and 
Sy) do not tend to the external Alfven frequency when d/R —> 2. 
The reason for this difference is probably that the A x and S y 
modes are strongly damped when d/R —> 2, and this fact has 
some impact on the real part of the frequency as well. 


On the other hand, Figure[2jb) shows that the high-frequency 
modes (A x and S y ) are more efficiently damped by resonant 
absorption than the low-frequency modes (S T and A v ). This 
result agrees with that of iRobertson & Rudermanl (201 ll) for 
large separations. However, the behavior of |o»i| /ojr obtained 
here for small separations is dramatically different from that 
of IRobertson & Rudermanl (201 ll) . They found that the oscil¬ 
lations become undamped in the limit d/R —> 2, while we 
find that the modes remain damped. Although the analysis of 
IRobertson & Rudermanl (201 ll) is mathematically correct, we 
find no physical reason for which the oscillations should be 
undamped when the loops are close to each other. Our com¬ 
putations show that the damping of the low-frequency modes 
(S x and A y ) is roughly independent on d/R, whereas the damp¬ 
ing of the high-frequency modes (A x and S y ) gets stronger as 
the se parat ion between loops is reduced. As pointed out by 
IRobertson & Rudermanl ( 2011 ) and iGiisen & Van Doorsselaerd 
(2014), the physical significance of the results obtained with by- 
cilindrical coordinates should be treated with caution when the 
separation between the tubes is small. 


Panels (c) and (d) of Figure[2]show the same results as panels 
(a) and (b), respectively, but now log 10 (l/R - 2) is used in the 
horizontal axes. These additional graphs are included to show 
in more detail the behavior of the solutions obtained with the 
T-matrix method for small separations between tubes. 

We have overplotted in Figure[2]the analytic approximations 
of wr and |o>i| /w R given in Equations (128b and ( 13 3k These ap¬ 
proximations were derived in the limit d/R^> 1 and reasonably 
agree with the numerical solutions when d/R > 3. As expected, 
the approximations do not work well for small separations. The 
analytic approximations were derived considering the contribu¬ 
tions from m — +1 alone, but the contribution of high m ’s to the 
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Fig. 2. Numerical results in the case of two identical coronal loops, (a) Dependence of cu R /w Ae on d/R, where tu A e = fc z v AjC is the external Alfven 
frequency, (b) Dependence of ]aii| /tu R on d/R. The meaning of the various lines is indicated within the panels. The dashed lines correspond to the 
analytic approximations in the limit d/R » 1 (Equations <28t and d33H . We have used f = 5, L/R = 100, and l/R = 0.2. Panels (c) and (d) show 
the same results as panels (a) and (b), respectively, but as function of log 10 ( d/R - 2). 




Fig. 3. Numerical results in the case of two identical coronal loops, (a) Dependence of w R /w A>e on L/R, where ru A e = fe,v Ae is the external Alfven 
frequency, (b) Dependence of |wi| /tu R on L/R. The meaning of the various lines is indicated within the figure. We have used £ = 5, l/R = 0.2, and 
d/R = 2.5. We note that the horizontal axes of both panels are in logarithmic scale. 


full solution is important for small separation between loops (see 
iLuna et aPl2009l) . 

Figure [3 displays the dependence of the solutions on L/R. 
This figure is included to show that the almost degenerate cou¬ 
ples S X -A y and A x -S y split into four different solutions for 
small values of L/R beyond the TT regime. We point out, how¬ 


ever, that the impact of the value of L/R on the solutions is 
not relevant when realistic values of this parameter are consid¬ 
ered. The TT limit used by iRobertson & Rudermanl (1201 ill and 
iGiisen & Van Doorsselaere! (i2014i) is therefore adequate. 

Now we plot in Figure [4] the dependence of wr and |tui[ /o»r 
on the nonuniform layer thickness, l/R. We consider a small 
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separation, namely d/R = 2.5, and the remaining parameters 
are the same as in Figure [2] Consistently, when 1/R = 0 the 
modes are undamped. The real part of the frequency of the 
low-frequency modes is almost independent of 1/R, while their 
|tui| / wr is roughly linear with l/R. Conversely, the real part of 
the frequency of the low-frequency modes decreases when l/R 
increases, and their |a>i| /wr is only linear with 1/R for small val¬ 
ues of this parameter. As discussed before, the behavior of the 
low-frequency modes (S x and A y ) is similar to that of the kink 
mode of an isolated loop. However, the high-frequency modes 
(A x and S y ) seem to be more affected by the interaction be¬ 
tween loops and show a somewhat different behavior when l/R 
increases. We note that because of the TB approximation we are 
restricted to consider small values of 1/R. 

It is useful to relate the present results with those of 
lArregui et al , (120071.120081) . who studied the damping of trans¬ 
verse oscillations of two nonuniform slabs. They found that the 
ratio IwiI/wr corresponding to the symmetric kink mode of 
the two slabs is weakly dependent of the separation between 
the two slabs (see lArregui et al.1 12007. their Figure 4). In turn, 
Arregui et al. (2008) found that the anti-symmetric kink mode of 
the two slabs is more efficiently damped than the symmetric kink 
mode (see their Figure 6). The symmetric and anti-symmetric 
kink modes of two slabs would be equivalent to the S x and A x 
modes of two cylinders. Thus, our results in cylindrical geometry 
are consistent with previous findings in Cartesian geometry. 

It is also c onvenient to consider the physical arguments of 
lArregui et al.l (120071 120081) to explain why the high-frequency 
modes damp more efficien tly than the low-frequency modes. 
lArregui et al . ( 1200 7, 12008 1 related the efficiency of the damp¬ 
ing with the magnitude of the total pressure perturbation within 
the resonant layers. According to lAndries et al.l (l2000l> . the ef¬ 
ficiency of the resonant coupling between the global transverse 
mode and the Alfven continuum modes is proportional to the to¬ 
tal pressure perturbation squared. ISolerl(|201 (>} ) plotted the square 
of the total pressure perturbation corresponding to the S x and 
A x modes (see his Figures 9.7 and 9.8) and found that, when 
the quantities are normalized, the perturbation of the A x mode 
reaches a larger value in the resonant layers than that of the S x 
solution. This result qualitatively explains why resonant damp¬ 
ing is more efficient for the A x mode than for the S x mode. 
Equivalently, a similar reasoning help us understand the different 
attenuation of the S y and A y modes. Nevertheless, a more robust 
study of the process of resonant absorption in two-dimensional 
configurations would be needed for a complete unders tanding of 
the different damping rates (see lRussett & Wrighlll20ldb . 

4.2.2. Non-identical loops 

Here we consider two loops with different properties and com¬ 
pare our results to those of iGiisen & Van Doorsselaere d2014i) . 
We use subscripts 1 and 2 to refer to the two different loops. For 
simplicity, we take R\ - R: = R and L/R = 100 in all following 
computations. 

First, we assume the same density contrast in the two loops, 
namely £i — £2 = 5, and vary I 2 /R while l\/R is kept fixed 
to 1\/R = 0.2. These resuits are shown in Figure [5] and can 
be compared to those already displayed in Figure [4] in the case 
of 1 1 = Z 2 . Importantly, we find that the collective oscillations 
remain damped when I 2 /R = 0. Although this is a rather par¬ 
ticular situation, the results have interesting implications. When 
L/R = 0 resonant absorption only occurs in the boundary layer 
of loop #1. However, this is enough for the collective oscilla¬ 
tions of the two loops to be efficiently damped. We note that in 


the study bv lGiisen & Van Doorsselaere] (i2014l ) the thicknesses 
of the nonuniform layers are linked to the coordinate system. 

Now we consider the case of two loops with different den¬ 
sity contrasts. We fix £2 = 3 and compute the solutions as func¬ 
tions of C\. These results are displayed in Figure [6j where the 
remaining parameters used in the computations are specified in 
the caption. For consistency, we keep the same notation as be¬ 
fore to denote the various modes according to the ordering of 
their frequencies, although they do not represent truly collective 
oscillations if £1 4 Ci ( Luna et al . 20091). Figure [6] can be com¬ 
pared to Figure 6 of Giisen & Van Doorsselaere ( 2014 ). To make 
aproper comparison, we note that IGiisen & Van Doorsselaerd 
(120141) plotted the signed damping ratio, while here we plot the 
absolute value. 

Concerning the real part of the frequency (Figure [6f a)), we 
find the same results as Luna et alJ (l2009t) . When £1 < £ 2 , the 
high-frequency modes are associated to loop #1 alone, whereas 
the low-frequency modes represent individual oscillations of 
loop #2. The opposite happens when £1 > £ 2 . Conversely, when 
£1 ss £2 the four modes approach and interact in the form of an 
‘avoided crossing’. Only in that case the modes represent truly 
collective oscillations (iLuna et alj|2009l) . 

Figure [6f a) also shows that the high-frequency modes are 
always within the Alfven continua of the two loops, i.e., the fre¬ 
quencies of the A x and S y modes are always larger than waj 
and co a .2 and smaller than oj Ac . This is true except in the limit 
£1 —> 1, where the frequencies of the A x and S y modes tend to 
tJA.e- On the contrary, the low-frequency S x and A y modes are 
below the Alfven continuum of loop #1 when £1 < 2 and below 
the Alfven continuum of loop #2 when £1 > 5. This result may 
have implications for the damping by resonant absorption. 

Figure [6fb) displays the damping ratio of the modes as a 
function of <fi ■ The result for the high-frequency modes can be 
understood as follows. When £1 —> 1 the frequencies of the 
high-frequency modes tend to a>A, e and, as a consequence, these 
modes become undamped in that limit. When 1 < £1 < ( 2 , the 
damping ratio increases when £1 increases. When 1 < £1 < £2 the 
high-frequency modes represent individual oscillations of loop 
#1. Then, the damping ratio reaches a maximum when £1 ~ £ 2 . 
When £1 > £2 the damping ratio saturates to a constant value 
because the high-frequency modes represent now individual os¬ 
cillations of loop #2, and the value of i /2 is fixed in the computa¬ 
tions. The behavior of the damping of the high-frequency modes 
agrees with that plotted bv lGiisen & Van Doorsselaera (2 0141) in 
their Figure 6. 

The overall behavior of the damping ratio of the low- 
frequency _modes displayed in Figure |6jb) also agrees with 
IGiisen & Van Doorsselaerd ( 2014 ). The damping ratio of the 
low-frequency modes is roughly constant when £1 < £2 and in¬ 
creases when £1 > £ 2 . The fact that the low-frequency modes are 
below the Alfven continuum of loop #1 when £1 < 2 and below 
the Alfven continuum of loop #2 when £1 > 5 have no important 
impact on the damping. Again, these results can be understood 
by considering that the low-frequency modes are associated to 
loop #2 when £1 < £ 2 , while they are associated to loop #1 when 
£1 > £2. 

The low-frequency modes computed here do not show the 
pronounced mini mum o f the damping rate seen in the solu¬ 
tion plotted by IGiisen & Van Doorsselaere (20141) when £1 ~ 
£ 2 . There are several effects that may explain this difference. 
The most_ obvious_one is_the different geometry considered 
in IGiisen & Van Doorsselaerd d2014l) and here. Another pos¬ 
sible explanation is that the density profile in the nonuni¬ 
form layers used by IGiisen & Van Doorsselaerd ( 2014 ) is dif- 
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Fig. 4. Numerical results in the case of two identical coronal loops, (a) Dependence of a> R /a) Ae on l/R , where w A _ e = Lv Ae is the external Alfven 
frequency, (b) Dependence of |tui| /wr on l/R. The meaning of the various lines is indicated within the figure. We have used £ = 5, L/R = 100, 
and d/R = 2.5. 



Fig. 5. Numerical results in the case of two non-identical coronal loops withfft = R 2 = R. (a) Dependence of a>R/a>A,e on L/R, where = k z v A , e 

is the external Alfven frequency, (b) Dependence of |tni| /wr on L/R. The meaning of the various lines is indicated within the figure. We have used 
ft = ft = 5, L/R = 100, h/R = 0.2 and d/R = 2.5. 




< i <, 


Fig. 6. Numerical results in the case of two non-identical coronal loops with R t = R 2 = R. (a) Dependence of a) R /a> Afi on ft, where u A ,c = k z v A,c 
is the external Alfven frequency, (b) Dependence of |tni| /% on ft. The dotted lines correspond to the Alfven frequencies of the two loops, and the 
meaning of the remaining lines is indicated within the figure. We have used ft = 3, L/R = 100, IJR = L/R = 0.2 and d/R = 3. 


ferent from that used here. A linear density profile is used in 
iGiisen & Van Doorsselaerej (l2014i) . so that the derivative of den¬ 


sity at the resonance position is independent of the frequency of 
the mode. Here we use a sinusoidal profile and take into account 
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that the position of the resonance and the value of the derivative 
of density at the resonance position are functions of the mode 
frequency. 

In Figure [b^b), the damping rate of the low-frequency modes 
shows a small bump around £) ~ 2, when the real part of the fre¬ 
quency approximately crosses the internal Alfven frequency of 
loop #1. The reason for this bump is that the the low-frequency 
modes intersect with and ‘avoid cross’ the fluting modes that 
cluster toward the int ernal Alfven frequency. This bump is ab¬ 
sent from Figure 6 of iGiisen & Van Doorsselaerel (120141) prob¬ 
ably because coupling between kink and fluting modes is not 
described in the TT approximation. 


5. CONCLUDING REMARKS 


In this paper we have extended the analytic T-matrix theory 
of scattering of iLuna et alJ ( 2009 1 12010 1 to investigate reso¬ 
nantly damped oscillations of an arbitrary configuration of par¬ 
allel cylindrical coronal loops. After presenting the general the¬ 
ory, we have performed a specific application in the case of two 
loops. This work is partially based on unpublished results in¬ 
cluded in I Soler ( 2010 ). where collective damped oscillations of 
prominence threads were studied. 

We have_compared our results to those of the pap ers by 
IRobertson & Rudermani d201 ll) and IGiisen & Van Doorsselaere! 
(120141) . They investigated the damping of collective oscilla¬ 
tions of two loops in the TT approximation and used a method 
based on bicylindrical coordinates. In general, the result^ of 
IRobertson & Rudermani ( 2011 ) and IGiisen & Van Doorsselaerj 
( 2014 ) are in good agreement with the present results, specially 
when the separation between loops is large. However, when 
the separation between the loops is small, i.e., for separations 
of few radii, the results of those previous works show impor¬ 
tant differences_compared to the present findin gs. F or instance, 
IRobertson & Rudermani (1 2011 ) and IGiisen & Van Doorsselaerj 
( 20141 ) obtained that by decreasing the distance between loops, 
the efficiency of resonant damping is reduced. In their compu¬ 
tations, both low- and high-frequency modes become undamped 
when the loops are in contact. However, this result lacks of a 
physical expla nation and contradicts pre vious findings in Carte¬ 
sian geometry dArregui et al.ll2007L 12008 ). In our computations, 
we find that the damping of the high-frequency modes gets 
stronger by decreasing the separation between loops, while the 
damping of the low-frequency modes is roughly independent of 
the separation. Our solutions do not become undamped when 
the two tubes are in contact. Thus , the results obtained here are 
consistent with previous results bv lArregui et al. ( 2007 1 12008 ) of 
collective oscillations of two slabs. 

Although the mathematical analysis of 
IRobertson & Rudermani (12011!) and IGiisen & Van Doorsselaerel 
(12014 is flawless, their results by may be affected by un¬ 
avoidable geometrical problems related to the bicylindrical 
coordinates when the loops are close to each other. In bicylin¬ 
drical coordinates the shapes of nonuniform boundary layers are 
not symmetric and change when the separation between tubes 
decreases. The nonuniform layers get thicker in the outer parts of 
the tubes and thinner in the inner parts. As already mentioned by 
IRobertson & Rudermani ( 2011 1 and IGiisen & Van Doorsselaerel 
(12014 . these geometrical limitations may lead to unphysi¬ 
cal results for small separations. The T-matrix method used 
here is not constrained by the geometrical problems of the 
bicylindrical coordinates. Therefore, we may conclude that the 
results given here are more generally applicable than those of 


IRobertson & Rudermani ( 2011 ) and IGiisen & Van Doorsselaerel 
(12014 when the loops are close to each other. 

Because of the TB approximation we were restricted to con¬ 
sider small values of l/R, i.e., l/R <sc 1. The effect of thick 
nonuniform layers could be included into the T-matrix formal¬ 
ism with the method of Frobenius used by ISoler et ali (1201 3l) in 
the case of an isolated loop. This would substantially increase the 
mathematical complexity of the problem but, on the other hand, 
it would provide a more accurate description of the damping 
of larg ely nonuniform loops. It has been shown by ISoler et al.l 
( 2014 ) that the error in the damping rate associated to the use 
of the TB approximation can be important when the loops are 
largely non-uniform. Apart from a numerical factor, the damp¬ 
ing rate in the TB approximation is independent of the specific 
density profile considered within the nonuniform boundary, but 
the density profile can have a more important impact when the 
non-uniform layers are thick. In addition, the real part of the fre¬ 
quency depends on l/R beyond the limit l/R <s 1. The effect 
of thick nonuniform boundaries on collective loop oscillations 
could be explored in the future. 

The method given here to compute resonantly damped col¬ 
lective oscillations can have multiple applications in the future. 
For instance, damped oscillations of a coronal arcade could be 
studied by modeling the arcade as a long line of parallel loops. 
Another interesting application is the investigation of oscilla¬ 
tions of loops formed by many strands. As shown bv lLuna et al.l 
( 2010 ), the global oscillation of the whole loop would be de¬ 
termined by the interaction of the oscillations of the individual 
strands. The resonant absorption process working in the individ¬ 
ual strands would affect the damping on the global loop motion 
(see iTerradas et al]|2008l) . In principle, the presence of multiple 
resonances in the system may cause the transverse oscillations 
of a multi-stranded loop to damp more quickly than the oscilla¬ 
tions of an equivalent monolithic loop. This idea could be con¬ 
firmed using the T-matrix method. Also, as shown bv ISoler et al.l 
( 2009 ), the effects of gas pressure and mass flow along the loops 
can easily be included in the T-matrix formalism, thus extending 
the applicability of the method. 
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